壁面函数在激波诱导分离流动中的应用

以数值模拟激波-附面层干扰引起的流动分离问题为研究背景,发展了基于有限体积方法的雷诺平均Navier—Stokes（RANS）方程的流场数值模拟方法。利用壁面函数模型得到壁面剪切应力,通过修正壁面粘性通量,构造了一种新的湍流边界处理方法,并将其耦合到RANS方程和SSTk-ω湍流模型的数值求解中;同时,针对激波诱导引起的附面层流动分离问题,提出一种附面层网格加密技术,能够自适应加密分离区内附面层网格,使得在流动分离区域也能够使用壁面函数模型。数值算例表明,壁面函数模型能够降低数值模拟结果对网格的依赖性;同时也验证了壁面函数耦合附面层网格自适应方法,在处理激波诱导引起的附面层流动分离问题时的有效性和准确性。     

三维烧蚀内部热响应数值计算研究

对热解型碳化复合材料三维烧蚀内部热响应数值计算关键技术进行了研究。采用碳化层—热解面—原始材料层模型,将热解气体与碳化层之间的对流换热处理为源项,通过有限元法建立移动边界条件下温度场求解方程组,采用Gauss Seidel迭代法计算热解气体质量流量和温度场。同时,研究和分析了三维烧蚀移动边界处理方法以及动网格生成方法。由于每个时间步都需要网格重划,烧蚀热防护数值计算对存储效率和计算效率要求较高,本文研究了内部热响应计算中影响存储效率和计算效率的主要因素,并提出了相应的压缩存储方案和求解方案。计算结果表明,移动边界处理方法准确合理;存储方案的存储效率较高;保持刚度矩阵和形函数矩阵正定对称性可以加快温度场计算的收敛速度。     

潜载导弹近筒口点火数值仿真

针对无横流影响的潜载导弹近筒口点火瞬态流场变化情况,建立多相流数学模型,考虑空穴模型,基于动网格技术,实现了潜载导弹弹射出筒、近筒口点火方式发射过程尾空泡的演变过程以及弹体、固壁的压强变化。仿真结果表明,筒口效应导致固壁各处的压强阶跃;近筒口点火后,尾空泡不断的颈缩、膨胀导致压强脉动,但不会产生断裂和回击现象。在筒口气团闭合后点火对筒口产生更大的压强脉动,有可能对固壁产生破坏。     

基于模糊理论的加速度计标度因数稳定性分析

运用模糊理论描述了影响加速度计标度因数稳定性诸多因素的不确定性问题,针对其中的磁钢磁感应强度和石英摆片热膨胀系数两类稳定性影响因素,给出了其三角模糊数区间,并通过该型加速度计的磁-结构耦合ANSYS分析求出了在模糊参数不同组合下的仿真输出,根据仿真输出结果获得了标度因数模糊稳定区间,给出了加速度计标度因数稳定性分析流程及实例.分析说明了基于模糊理论的稳定性分析方法可以为加速度计加速稳定试验截止条件提供依据,对提高加速度计标度因数长期稳定性水平具有重要意义.     

一种基于三角形重投影的虚拟视角合成

相比较传统的视角插值合成方法只能在基线上合成虚拟图像,提出的一种基于三角形重投影的虚拟视角合成方法,可以在已知少数几个匹配点的情况下,将场景内容自动迭代划分为一致和非一致三角形区域,通过不同策略的三角形重投影实现更加自由的视点的虚拟图像合成.划分过程中引入直线段约束可以解决直线纹理在三角形接缝处的畸变问题,保证最后的合成图像中场景表示在此处的连续性.利用设计的可信度函数,沿一定搜索路径可找到更多的新匹配点,从而带来更多的一致三角形区域.将此方法应用到虚拟彩色图像的合成,实验结果证明了此方法可以在更大视点范围内合成逼近于真实图的虚拟图像.     

非结构变形网格和离散几何守恒律

数值模拟流固耦合问题或多体分离问题的非定常流动时,常采用基于任意拉格朗日-欧拉（ALE）方程的有限体积法,涉及到变形网格和离散几何守恒律。在对变形网格算法进行综述时,按照构造思想分为物理比拟、椭圆光顺、插值、运动子网格（MSA）及其混合法共5类,分别介绍了基本原理、研究现状和适用范围,通过算例比较表明,径向基函数（RBFs）和运动子网格相结合的混合方法既有很好的变形能力,也有较高的计算效率,值得进一步发展和推广。在介绍了离散几何守恒律（DGCL）概念之后,采用二维几何模型进行分析,指出其机理是离散过程中体积增量与网格面元扫过的体积不相等造成的,把目前国内外应用的算法分为面积修正法、给定速度的面积修正法、速度修正法和体积修正法共4类,对其应用范围和存在的问题进行讨论,认为提出的体积修正算法既可以保证流固界面条件,也可以用于时间多层格式。     

基于非结构网格的离子推力器羽流数值模拟方法研究

为适应复杂的航天器外形,文章采用基于非结构网格的数值模拟方法研究离子推力器真空羽流的基本特征。首先,介绍了非结构网格划分的方法,并分析了建模区域生成各种拓扑关系的特点;然后,提出了粒子穿越网格的具体算法,包括有利于减少计算量的结合矩阵重标号的带宽压缩算法;最后,开展了模拟结果和试验结果的对比研究,验证了以上数值模拟方法的可行性。     

固体姿轨控针栓推力器内流场数值仿真

为研究针栓推力器内.流场的特征和规律,采用二维轴对称模型,开展了针栓运动过程的非稳态及稳态流场仿真,分析了不同开度下的流场结构、针栓所受气动力及流场对针栓几何参数的敏感性.结果表明:由于针栓的干扰作用,流场结构呈现出复杂激波系、顶端流动分离等现象;气动力系数与针栓头部锥角大小呈正相关关系,与集气室压强无关;节流面形状对气流参数分布、激波强度及位置、气动力等表征的流动模态影响很大.     

用于气动声学计算的非均匀网格紧致差分格式

为克服传统紧致差分格式在数值求解非均匀网格问题时产生的寄生波,构造了一种新的高精度紧致差分格式.通过泰勒展开分析方法,详细给出了格式系数的通用形式;利用傅里叶分析方法,分析了数值耗散、色散误差.以3对角6阶精度紧致差分格式求解均匀扰动网格问题为例,计算表明:色散值和耗散值随扰动因子的增加而更加趋近于精确值;当扰动因子大于0.213时,格式不稳定,当扰动因子小于等于0.213时,格式渐近稳定;对一维对流波和二维波传播的模拟计算所得数值解与精确解吻合,体现了该格式在求解非均匀网格问题时的优越性.      

基于非结构嵌套网格方法的旋翼地面效应数值模拟

 建立了一个基于非结构嵌套网格的流场求解器,用来精确模拟复杂的旋翼近地流场,为更好地分析地面效应(IGE)对旋翼气动特性的影响提供一套计算方法。在该求解器中,控制方程采用惯性坐标系下的非定常N-S方程,空间方向上采用二阶迎风格式,并用背景网格的一个面模拟地面作用,以方便地面边界条件的处理。应用所建立的模型,首先针对有实验结果可供对比的旋翼无地面效应（OGE）流场进行了数值模拟,以验证计算方法。然后着重计算了IGE下流场中的桨尖涡空间位置和旋翼拉力增益,并对旋翼在小速度前飞状态下的近地流场及地面涡形成过程进行了模拟和分析。在此基础上,得出了一些有意义的结论。